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Introduction 

A number of global atmospheric general circulation models (GCM's) 
have been developed during tiie past 20 years with the objective of 
accurately simulating the large-scale dynamics and physics of the 
atmosphere. In some CCM experiments the models have been started from 
a hypothetical barotropic state of rest and allowed to "spin up" to a 
climatology comparable to that of the present terrestrial atmosphere. 

In others they have been initialized with real, or realistic, data and 
allowed to generate a simulated forecast meteorological history. In 
either type of experiment, the validity of the model is tested by 
comparing the model -generated climatology or forecast with the real 
atmosphere. These comparisons have been carried out through diagnostics, 
sucli as gross energetics, transports, momenta, and hydrologies, as well 
as vertical -meridional cross sections of various zonal mean atmospheric 
properties. Model and nature have also been compared in terms of 
horizontal synoptic fields, using traditional measures of forecast skill, 
such as those commonly employed for the verification of prognostic weather 
maps . 
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No method of verifying prognostic maps (whether applied to explicit 
daily forecasts, monthly mean forecasts, or simulations of climatology) 
is completely satisfactory. Subjective comparisons of synoptic patterns 
are obviously inadequate, and quantitative scalar measures of agreement 
(e.g., root-mean-square errors, correlation coefficients, and gradient 
"skill scores"), while they may be objective, convey little information 
about pattern agreement, except in a gross relative sense. Additional 
information may be provided by a comparison of the spectral components 
of the patterns, as represented, for example, by the expansion coefficients 
of a series of surface spherical harmonics. 

Spherical harmonic analysis has been applied in the past to many 
kinds of geophysical problems (e.g.. Chapman. and Bartels, 1962), and is 
now even used operationally for global weather analysis at the National 
Meteorological Center (NMC) . However, it has apparently not been widely 
adopted as a method of synoptic pattern verification. (Bor an example, 
see Leith (1974).) In this study, surface spherical harmonics are 
used to analyze the horizontal fields of various quantities generated 
by one global GCM - tiie C1SS "climate model" (Hansen et al., 1979) - and 
to compare the model results with nature. 

Forecast experiments have been carried out with the climate model 
by initializing it with global NMC "data" (actually derived from NMC 
operational analyses) for the first day of a month (at 00 GMT) , running 
the model to the end of the simulated month, and computing the monthly 
mean fields of various predicted quantities. These monthly mean forecast 
fields arc then compared with the corresponding observed fields for the 
same month. The quantities selected for analysis arc the sea-level 
pressure (SLP) in millibars (mb), the 850 mb temperature (T8) in Kelvin 
degrees (K) , and the 500 mb geopotential height (Z5) in meters (m) . 
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Monthly mean (observed) climatologi . al fields of these prcdictands, as 
well as model -generated "climatologies", are also used in the evaluation 
program . 

Any global horizontal field (forecast, observed, or climatological) 
can be expanded into a finite series of spherical harmonics. (For further 
details on applications of spherical harmonics, sec, e.g., Chapman and 
Bartels (1962), Spar (1950), Belousov (1962), Merilees (1973), and 
Blackmon (1976).) The fields may be represented by tables of the 
expansion coefficients, or by tables of the magnitudes and phases of the 
spectral components, which may also be shown, for selected wave numbers, 
in the form of vector diagrams ("harmonic dials"). Verification of the 
model output can then be expressed in terms of the relative magnitudes 
of the dominant harmonics, with phase angles also considered. (The 
practical advantages and limitations of spherical harmonics, as compared 
with Fourier series and Chebyshev polynom ials * for the solution of 
atmospheric problems on the sphere arc discussed by Boyd (1978).) 


Spherical harmonic expansion 

A horizontal field, Q, may be expanded as a function of latitude, , 
and longitude, X , into a series of surface spherical harmonics, 

• t t (C... Y“ . $,. Y “ ) CO 

where Yn,in are the even (e ) and odd (o) normalized spherical harmonic 


functions (see, e.g., Blackmon, 1976), 
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"zonal harmonics"), Pn ,m(sin $> ) are the associated Legendre polynomials 
("spherical surface harmonics"), m is the zonal (longitudinal) wave 
number, n-m indicates the number of modal parallel'; (twice the meridional 
wave number) from pole to pole, N is the truncation degree of the 
series, and Cn,m and Sn,m arc the normalized expansion coefficients of the 


series, to be determined from the data, Q, by 

Is'i-TfM 
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The orthonormali ty property guarantees that the area-weighted mean 
value of (Yn,m* Yn' ,m' ) over the sphere, for either (c) or (o) , equals 
zero if n + n' or m ^ m' and (4“rr) -1 if n = n' and m = m' . This is the 
basis for (3) and the calculation of the coefficients. Normalization 
makes it possible to compare coefficients of different order, m, and 


degree, n, in the same series or among different series. 

? 2 2 

The magnitude, An, in, of any harmonic is given by ATi,m = Cn,m +Sn, m, 
the phase angle , P by * tan"* Sn,m/Cn,m, as indicated in the 

• »Vl 

harmonic dial in figure 1. 
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Various numerical integration schemes for the calculation of (3) 
were tested by checking both the orthonormality properties and the 
errors of reproduction of given fields. Also, tests were carried out 
with different degrees of truncation, the maximum being N = 18. 

(Alternative numerical methods for computing associated Legendre functions 
and expansion coefficents are discussed by Mcrilecs, 1973.) The method 
finally adopted was found to he satisfactorily accurate when applied to 
data fields on either an S° x 10° or 4° x 5° latitude-longitude grid. 

While tables of normalized associated Legendre polynomials arc 
available (c.g., Belousov, 1962), they generally do not correspond to 
the latitude interval of the HISS grid. Therefore, for the present 
study, the polynomials were calculated for tne appropriate latitudes 
by means of Rodrigues' formula (Hobson, 1955), and stored on disk for 
use as required. A comparison of computed and tabulated polynomials 
for the same latitudes up to degree 18 gave perfect agreement to the 
10 -th decimal place. 

The coefficients, Cn,m and Sn,m arc calculated from (2) and (3) 
in steps, integrating first over <|>» then over A, using a combination of 
the Simpson, "three -eighths" , and trapezoidal rules, with a smoothing 
operator applied at v> > 7 , for the ^ "integration, and Filon's rule 
(Davis and Kabinowitz, 1967) for the A - integrals. 

The area -weighted mean value of Q over the sphere is given by 
(4lf ) */ 2 Co,o. In general, Cn,o, the coefficient of the one-dimensional 
spherical harmonic of degree n ( a function witli n nodal parallels), 
represents the amplitude of a zonally symmetric function of latitude 
only. For example, Y 1>0 = sin<|?, Y, 0 = (3.2 T)* (3 sin 2 ^-l), 

etc. Thus, the difference between the values of Q at the North and 
South Foies is represented by approximately ^l»o a, 'd the difference 
between the polar and equatorial mean values of Q is represented 


- 5 - 


by approximately C-, 0 . Polynomials of even degree are symmetric 
about the Liquator, while those of odd degree arc ant i -symmetric . The 
coefficient C^q, for example, represents a component that is symmetric 
about the Equator, with a minimum (or maximum) at the Equator and three 
alternating bands of high and low values in each hemisphere (as in the 
case of the planetary sea-level pressure distribution). 

The associated Legendre functions represent two - dimensional wave 
patterns. Tlius, for example, ^ represents the checkerboard pattern 
illustrated in figure 2, with 2 nodal parallels, 2 longitudinal waves, 
and alternating high and low values. 

Some Illustrative Examples 

Synoptic monthly mean maps of SLP, T8, and Z5 (observed, climatological, 
and forecast) are to be expanded in finite series of spherical harmonic 
functions up to N=18. The observed monthly mean fields (0) are derived 
from NMC operational analytic data obtained from the National Center 
for Atmospheric Research (NCAR) and interpolated to the G1SS climate 
model grids. Monthly mean forecasts (F) are computed with the climate 
model, starting from OOGMT initial conditions on the first day of the 
month. The actual climatological (C) fields for each calendar month 
were provided by NCAR. (A set of monthly mean "model climatologies" (M) 
is also being generated by running the medium mesh S°xlO° , climate model 
for 5 simulated years and averaging each month's output. Evaluation of 
the model will be carried out in terms of both its climatology error, 

M-C, and its "anomaly error", (F-M) - (0-C).) 
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To facilitate interpretation of the spherical harmonic expansions, 
the leading harmonics for each field are tabulated and ranked in 
descending order of the magnitules of the amplitude coefficients An,m. 

Also listed, for m > o, are the phase angles. The sign of An,o in a 
table is actually the sign of Cn,o. Synoptic maps may then be compared 
in terms of the coefficients and phases of the dominant components of 
each field. Area-weighted global mean values (from Co,o). are also 
given for each field. 

Table 1 shows some sample results for the December climatology (C) , 
the December 1976 observed (0) fields, and a forecast for December 1976 
(F) made with a 7-layer medium mesh (8°xl0°) model (MX444M7). Only the 
8 largest harmonics are listed for each field, together with the 
global mean values. 

The data listed in Table 1 indicate that reasonable global mean 
values are "forecast" by the model, but with a slight cold bias. The 
same dominant harmonic (2,0 for T8 and Z5 and 1,0 for SLP) appears in 
the fields of C, 0, and F. However, the model (F) overestimates the 
polar-equatorial difference (2,0) in T8 (41.5 vs. 37.8) and underestimates 
it in ZS (800 vs. 855), while the interhemispheric difference (1,0) in 
SLP is grossly underestimated (11.9 vs. 22.3). 

The differences between the 8 leading harmonics for C and for 0 
indicate the character of the monthly anomaly. For example, in the 
expansion of T8, 6 of the 8 leading climatology harmonics are represented 
among the 8 leading observed harmonics. However, two outstanding 
observed anomalies are the (6,0) and (4,1) components, which do not 
appear among the first 8 climatology terms. One of these harmonics, 

(6,0), does appear among the leading forecast terms, albeit too large. 

On the other hand, the model generates some components (3,0 and 6,2) 
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Table 1. Global mean values and the 8 largest harmonics (An,m) for each field (T8, Z5, SLP) for 
December, climatology (C) , and December 1976 observed (0) and Forecast (F) , medium 
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which do not appear among cither the 8 leading obser’ ed or climatological 
harmonics . 

Similar results are indicated for Z5 and SLP. In the case of Z5, 
the anomaly (difference between 0 and C) is represented by the last 4 
harmonics in 0, none of which appear among the. 3 leading C coefficients. 
The forecast error (difference between F and 0) is represented by the 
forecast terms, (4,0), (5,0), (6,2) and (4,2), which do not appear in 
the first 8 observed harmonics, and the observed coefficients, (4,1), 
(7,4), (6,5), and (5,2), which are not found among the first 8 forecast 
terms. In SLP the anomaly is represented by (6,1) and (8,1), which 
are present in 0 but not in C, and by (4,0) and (5,0), whicli appear in 
C but not in 0. The forecast error is dominated by (4,0), which is not 
found in the observed list, by (6,0), which does not appear among the 
first 8 forecast components, and by (2,0), for which the forecast and 
observed harmonics are opposite in phase. The mismatch between the F 
and O harmonics is particularly evident in SLP. 


- 8 - 


References 


Belousov, S. L., 196.’: Tables of Normalized Associated Legendre 

Polynomials (translated by D.E. Brown). New York, Pergamon 
Press, 379 pp. 

Blackmon, M.L., 1976: A climatological spectral study of the 

500 mb geopotential height of the Northern Hemisphere. 

J .Atmos . Sci . , 33, 1607-1623. 

Boyd, J.P., 1978: The choice of spectral functions on a sphere 

for boundary and eigenvalue problems: a comparison of 

Chebyshev, Fourier and associated Legendre expansions. 

Mon. We a. Rev., 106, 1184-1191. 

Chapman, S. and J. Bartels, 1962: Geomagnetism . (Volume II. 

Corrected edition.) London, Oxford University Press, 1049 pp. 
(First edition, 1940) 

Davis, t'.J. and P. Rabinowitz, 1967: Numerical Integration . 

Waltham, Mass., Blaisdell Publishing Co., 230 pp . 

Hansen, J. and Collaborators, 1979: An efficient three-dimensional 

global model for climate studies. Goddard Institute for 
Space Studies (GISS) . (Unpublished) 

Hobson, L.W., 1955: Theory of Spherical and Ellipsoidal Harmonics . 

New York, Chelsea Publishing Co . , 500 pp. r— 

Leith, C.E., 1974: Spectral statistical -dynamical forecast 

experiment. International Symposium on Spectral Methods of 
Numerical Weather Prediction, Copenhagen, August 1974. 

GARP Rep. No. 7. 


Mcrilees , P . L . , 
of a series 
224-227. 


1973: An alternative scheme for the summation 

of spherical harmonics. J. Appl. Meteor, 12, 


Spar , J . , 1950 : 
J. Meteor., 


On the theory of annual pressure variations. 
7, 167-180. 


- 9 - 



